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Abstract 

The use of correlation between two open-charm mesons is suggested to give 
information about the nature of the medium created in heavy-ion collisions. In- 
sensitivity to the charm production rate is achieved by measuring normalized 
cumulant. The acollinearity of the D momenta in the transverse plane is a mea- 
sure of the medium effect. Its dependence on nuclear size or Et provides a 
signature for the formation of quark matter. 
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The conventional probes of quark matter in heavy-ion collisions, such as dileptons and 
J /if) suppression, have not provided conclusive evidence about the creation or absence of 
quark-gluon plasma []]]. A large part of the difficulties involve the ambiguities arising from 
competing processes and from uncertainties in the initial normalizations of some key quan- 
tities. An effective probe should be free of such ambiguities. In this paper we suggest the 
possibility that charm correlation may be such a probe. 

The correlation proposed is between open-charm mesons at nearly opposite directions. 
The charm quark is used both for tagging and for probing. Heavier quark can be used 
when appropriate; charm will be used as a generic term for heavy quark in the following 
discussion. Lighter partons and associated minijets are too copiously produced at RHIC 
and LHC [0 to be useful for our purpose of tagging and probing. The idea is based simply 
on the dual requirements that the signature should be independent of the production rate 
but sensitive to the medium through which the probe traverses. Appropriately normalized 
cumulant can satisfy the first requirement, while the transverse deviation from exact back- 
to-back correlation meets the second. 

Briefly stated, it is suggested that one searches for DD produced in the transverse plane at 
y = with their momenta nearly collinear, but opposite. The acollinearity in the transverse 
plane is the measure of interest. To enhance the effect, experimental cuts should be made 
on the magnitudes of the D-meson momenta so that they are nearly equal and not too large. 
We expect that the mean acollinearity is smaller if the medium is deconfined quark matter 
than if it is not. 

The proposed measure is similar in spirit to the acoplanarity of jets suggested by Appel 
0, but significantly different in substance. The major differences are: (a) the jet axes 
cannot be as precisely determined as the D-meson momenta, (b) at high energy too many 
jets are produced resulting in contamination and deterioration of the correlation signal, (c) 
we emphasize the difference between the propagation of a c quark through a deconfined 
medium and that of a D-meson through a confined medium, and (d) the phenomenology of 
DD correlation can reveal interesting physics even in kinematic regions where perturbative 
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QCD (pQCD) is unreliable. 

Let pi and p 2 be the momenta of the two detected charge-conjugate D mesons. The 
cumulant is 

c(Pi,p 2 ) = toipxiih) - Pl(Pl)Pl(P2), (1) 

where p n is the n-particle distribution function; it is the irreducible part of the two-particle 
correlation. In a heavy-ion collision the cumulant can in general be expressed in the form 

c(pi,P2) = I d ^ d J^ S(kx, k 2 )H{k 1 ,p x )H{k 2 ,p 2 ), (2) 

where S(k\, k 2 ) is the probability of producing two partons with momenta k± and k 2 , and 
H(ki,pi) is the hadronization function that connects the parton i at the point of creation 
to the hadron detected with momentum pi. Hereafter we adopt the convention of using the 
symbol k (p) for parton (hadron) momentum. It should be stressed that H is not simply the 
fragmentation function usually used for jet considerations because firstly the produced parton 
must traverse a dense medium and suffer momentum degradation before fragmentation, and 
secondly the hadronization process may be recombination [|], [5] instead of fragmentation. 
In fact, it has been shown that the data of open charm production in the forward region 
of hadronic collisions can be well described by recombination M, but badly by pQCD or 
fragmentation model M. 

In the domain where pQCD is reliable one can write S(ki, k 2 ) for AB collision as 

S(h, k 2 )=cj d 2 ad ^ b F A (k a )F B (k b )S 4 (k a + k b -h- k 2 ) \M(a + b^l + 2)\ 2 (3) 
J A;" k u b 

plus other terms of similar structure, if more than one hard subprocess are important. In (^) 
c is a numerical constant, FA(k a ) is the parton distribution in nucleus A, and M(a+b — > 1+2) 
is the amplitude of the hard subprocess involved. There is a great deal of physics contained 
in the determination of Fa,b, which depends on nucleon structure function at small x, gluon 
distribution, nuclear shadowing, initial-state radiation, preequilibrium and possibly thermal 
interactions, space-time evolution, etc. So much uncertainty is involved in the problem that 
the study of open-charm production has been suggested as a means to learn more about the 
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parton dynamics in the early phase of nuclear collision [KL |9|, i.e., the reverse of using Fa,b to 
predict measureable quantities. While that is certainly a worthwhile project to pursue, our 
proposal here is to circumvent all that complication and proceed with the use of S(ki,k 2 ) 
independent of the details about Fa,b- 

From S(ki,k 2 ) n °t only can two-particle inclusive distribution be determined as in 
the one-particle distribution can also be obtained as follows 

PM) = I ^^S{k 1: k 2 )H{k 2: p 2 ). (4) 

To free our signal from the uncertainties of the primordial parton dynamics, let us define the 
singly-normalized cumulant function 

C(pi,P2) = c(pi,p 2 )/pi(p 2 ). (5) 

It is clear from (0) and (f|) that C(pi,p 2 ) should be insensitive to the rate of charm produc- 
tion. 

Hereafter we shall regard hadron 2 (with momentum p 2 ) as the trigger particle, against 
which we study the properties of the probe particle 1. Of course, it is the relative momentum 
between the trigger and probe that is important, but conceptually it is efficient to identify 
(arbitrarily) one of the two D mesons as the trigger and define the axes such that the trigger 
momentum in every event is always aligned along a fixed direction, say — x axis, with y being 
the axis normal to the scattering plane containing the beams and the trigger. The aim is to 
study the momentum distribution of the other D meson in the neighborhood of the +x axis. 

Since the parton momenta k\ and k 2 can have large longitudinal imbalance due to unequal 
momenta, k a and k^, of the initial colliding partons, but they can have only limited total 
transverse momentum Kt = k\ T + k 2r due to the small intrinsic transverse momenta of 
the initial partons and to the initial-state radiation, we can avoid the complexity of the full 
structure of S(ki,k 2 ) if we restrict k\ and k 2 to only the near neighborhood of a common 
transverse plane. That is achieved by requiring that p\ and p 2 lie only in the transverse 
plane at y = 0. For brevity we shall refer to that plane as Tq. Since we expect the angular 
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differences between k{ and pi to be small, that requirement therefore forces k± and &2 to be 
very close to % also. 

The aim of this problem is to learn about the medium effect through H(ki,pi), on which 
the measureable C(pi,p2) depends. With k{ near 7o those partons do not participate in the 
longitudinal expansion of the system, which is another area of large uncertainties. But even 
in % there are several possible processes leading from the partons (ki) to the hadrons (pi), 
each involving a different hadronization function H{k, h p,j). So far we have not specified the 
kinds of partons carrying fcj. They can be high- momentum quarks of the u and d types, 
or lower-momentum gluons, all capable of fragmenting into the D mesons. Since the two 
fragmentation processes are independent, p\ and pi are not correlated even though k\ and 
^2 are. To narrow down the hadronization process we make use of the experimental freedom 
to require further that the magnitudes p\ and P2 are nearly equal within a narrow range. 
Moreover, that magnitude should not be too high, say, in the 2 — 5 GeV range. In so doing 
we can maximize the contribution from cc pair creation to the formation of DD. There are 
several stages of reasoning involved here, which we now describe. 

The possible parton types are q and c, where q denotes u, d, or g collectively. Let the flavor 
labeling not be encumbered by concerns about quark or antiquark differences. We postpone 
our consideration about the s quark until later. The possible hadronization processes are 
fragmentation (F) and recombination (R), for which k > p in F, but k < p in R. Thus for 
the production of D there are four possible processes: F(q — > D), F(c — > D), R(q — > D), and 
R(c — > D). At high collision energies there are so many hard subprocesses that there are 
enough transversly moving partons to make recombination competitive with fragmentation 
in the formation of D in T . That is not the case in pp or pA collisions. In fact, for any 
given momentum p of D in %, it is more favored to recombine two lower momenta partons 
to add up to p than to create a higher momentum parton which subsequently decays to p, 
since the probability of creating high kr partons falls off as a power k^ a , with a > 4. For 
single-particle inclusive distribution a comparison between the two hadronization processes 
in the production of particles in % has been studied quantitatively with the result that 
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for pt < 6 GeV/c and for large nuclei the rate of hadronization through recombination is 
at least an order of magnitude higher than through fragmentation. For a pair of correlated 
particles the R/F ratio of the rates would be squared. Thus we may ignore F(q — > D) and 
F(c — > D) in the following discussion. 

Since m c rs 5 m, for constituent quark masses, the momentum fractions x c and x q of 
c and g, respectively, in D are on the average very different, with x c ~ 5x g ; hence, the 
recombination function for c + q — > D is maximum when k c m bk q in the same direction 
p. . For p*! and p 2 nearly equal and opposite, the production of DD is therefore dominated 
by the creation first of cc pair with k\ ~ — &2 followed by recombination with low-momentum 
q quarks, rather than by the process where a created qq pair dictates the momenta of the D 
mesons. In short, R(c — > D) is more important than —►£)). In the following we shall 
focus on the process where S^fci, ^2) m (@) describes the hard production of a cc pair, and 
the two H functions represent R(c — > .D). 

The only part in the problem that has a firm theoretical footing is the amplitude M 
for hard scattering in (Q), which is calculable in pQCD. Even there, charm production 
with kr ~ 2 GeV/c is in the grey area of reliability. As mentioned earlier, the parton 
distribution F^(x a ) and F B (xb) are quite uncertain, but the normalized cumulant C(pi,p2) 
is insensitive to all of them, including M; furthermore, we emphasize the misalignment of 
Pi and f>2 , their magnitudes being selected by experimental cuts. Thus the signature we 
seek depends mainly on H{ki — > pi), which is sensitive to the medium that stands between 
the creation of cc and the detected hadrons DD. That is just what a good probe should 
be. It is unfortunate that H(ki — > pi) cannot at this point be calculated precisely in QCD, 
perturbative or otherwise. However, the discovery of unambiguous experimental signature 
is more important than having reliable theoretical calculations at this stage. On the basis 
of reasonable arguments we indicate below what that signature might look like. 

The four vectors k{ and pi (i — 1,2) are all very close to %. We consider below only 
their projections k iT and p iT on %. For brevity we shall omit the subscripts T, unless 
there is confusion. Let the angles among these four vectors be labeled as follows: 6' 1 (/c 1 ,p 1 ), 
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#2(^2, P2), 6*12(^1,-^2), and 4>(pi,— P2). For simplicity let the distributions in these angles 
be represented by Gaussians: exp(— 2 /A 2 ), exp(— #2/A 2 ), exp(— 6^/A 2 ^), and exp(— 2 /A 2 ), 
respectively. The distribution in is a Gaussian also because of the convolution theorem, 
since (0) implies 

c(pi,p 2 ) oc J d6 1 d6 2 exp I —± — — -^§1 , (6) 

where (9 12 = <fi — 61 — 9 2 - Thus if Ai, A 2 and Ai 2 are all small compared to ir, then 

A 2 = A 2 + A 2 + X\ 2 . (7) 

One can improve on these integrations with better knowledge on S{k\, k 2 ) and H(ki,Pi), but 
the result will not differ much from (^), which encapsules the property that the width in 
PfP2 is the rms sum of those in pyki , kvk>2 and p 2 -fc 2 . Since we require \p\\ ~ \p\\ = p, we 
may write C(pi,p 2 ) as C(p,(f>). Then we have 

C{pA) = C{p)e-^l x \ (8) 

Alternatively, one can define p\ = — p 2 + p t (note the small t), where p t is nearly normal to 
the trigger axis but still in %. Thus experimentally, the data on C(pi,p 2 ) can be presented 
as distributions in £ = p t /p for various values of the trigger momentum p. Empirical C(p, £) 
should be sharply peaked in £, though not necessarily Gaussian; the width can be represented 
by A in @ in the following discussion. 

In order to keep all the quantities in (0) small, p should not be too small; neither should it 
be too large so that the rate of producing two correlated D mesons does not become too low. 
The range 2 < p < 5 GeV/c appears to be reasonable. In that range the c and c quarks are 
beams in the dense medium with momenta k\ and /c 2 , both < p, insufficient to give any hope 
to the validity of pQCD in describing their passages through the medium. Nevertheless, we 
expect charm quarks in that momentum range to be sensitive to the medium effects and can 
provide us with useful information. 

We now consider the various contributions to (^|). Nonzero #i 2 means that K = k\ + fc 2 ^ 
0. Nonzero K is due to the intrinsic transverse momenta of the partons plus the recoil from 
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initial-state radiation before the hard subprocess that creates the cc pair. We expect \K\ to be 
small, and set its average value at 0.3 GeV/c for low fcj and allow it to be higher for higher ki. 
Approximating k\ = k 2 and denoting them collectively as k, we have cos 9 12 = 1 — K 2 /2k 2 , 
and 9 12 — K/k. Thus for k > 2 GeV/c, we have A12 < 0.15. More importantly, A12 is 
independent of the medium in which the created c quarks will traverse, and will therefore 
not affect our signature of the medium effect. 

Next, we consider the last two terms of (^j) due to the H functions in (||), which are 
the heart of the problem. The two H functions independently describe the hadronization 
processes from ki to pi, as they proceed along the paths l\ and l 2 , through the medium in 
essentially opposite directions in % from the point of creation of cc. Since the result we 
seek concerns the sum A^ + A| in (|?|), it is equivalent to that due to one c quark traversing 
the entire path L = 1% + l 2 from one end of the medium to the other in %, passing through 
the point of cc creation. The transverse expansion of the medium is not rapid. If R is 
the average radius of the system that is relevant to this study, then when averaged over all 
possible points of the cc creation in T and over all orientations of the path, the mean path 
length L of crossing a disc of radius R is L = 4R/3. Our consideration is thus reduced to the 
hadronization of a c quark with initial momentum k passing through a slab of the medium 
of thickness L and emerging as a D meson with momentum p. 

There are two scenarios to examine. One is that the medium consists of deconfmed quarks 
and gluons, while the other is of high-density hadrons with quarks confined. Let them be 
referred to as quark matter (QM) and hadron matter (HM), respectively. Of course, they 
represent two extreme cases, and other scenarios that stand between them are possible. If we 
know the nature of the signatures for the extreme cases, what happens in the intermediate 
cases can be estimated by interpolation. Thus for now it is sensible to consider just QM and 
HM. 

1. Quark Matter. Being deconfined, the medium cannot support the formation of any 
hadrons nor the existence of any color flux tubes in it. Thus the c quark that traverses the 
QM medium remains as a c quark. It may lose momentum and deviate from straightline path, 
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but the formation of D can occur only after passing through the medium by recombining 
with a q at the exit point. If pQCD were applicable, one could study in detail the effect 
of multiple scattering as in |[T[ [1^], and determine the degree of energy loss and k t gain. 
Hereafter we use k t to denote parton momentum transverse to k. However, for k as low as 
2 GeV/c, the reliability of pQCD is questionable. Qualitatively, one expects the radiative 
energy loss to be reduced for heavy quarks compared to light quarks. We assume that the 
loss of longitudinal momentum is roughly compensated by the gain in momentum due to 
recombination so that p ~ k. The more important aspect of the problem is the k t gain. If one 



regards the result of [TT| that takes the LPM effect into account as being valid, then k t due 
to gluon radiation is of the order of the color screening mass fi, independent of the number 
of multiple scatterings. The cumulative effect on k t due to elastic scatterings depends on 
whether the random-walk model is valid or the quantum coherence effect is important. In 
the former case the process is Markovian and k t would increase with L, while in the latter 
case it would not. The coherent LPM effect in the longitudinal component is non-Markovian. 
It has been known that the former is more relevant to conventional large-p^ processes, while 
the latter is for low-pr processes. Our H(c — * D) in the present problem belongs to the 
latter category. The following experimental facts support the latter. 

Exhaustive studies of h\A — > h 2 X inclusive reactions at high energies and low px have 
revealed that, for hi in the beam fragmentation region of h%, the pt distribution of hi is 
essentially independent of A |13|, [14]]. Because of the flavor dependence of the hi distribution, 
the processes can be well interpreted in the parton picture by considering the valence quark 
that is common in h\ and hi ||, |H], [L5[]. Take h\ = tt + and hi = K + to be specific. For K + 
with high xp, it is the high momentum u quark in ir + that leads to K + by recombination. 
The pt of K + reflects the &t of the u quark after traversing the target. The independence of 
(Pt) on A implies that there is no kj- broadening due to multiple-scattering effect on the u 



quark. Specifically, the data of |T^] indicate that a(pr = 0.5)/a(pT = 0.3) stays essentially 
uncharged at 0.5 and 0.48 for A = Cu and Pb, respectively, for p n + = 100 GeV/c and 
Pk+ = 80 GeV/c. Although it is a u quark traversing normal nuclei, whereas our problem 
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involves a c quark traversing dense QM, the independence of kt on the path length is likely 
to be due to a common origin: quantum coherence in low-k t processes is non-Markovian, so 
kt does not increase with L. The same lack of A dependence is found in |13[] for hi = p and 
hi = A. We emphasize that it is true only at low px- Significant A dependence at high pt is 
not excluded, such as in the production of massive dileptons in pA collisions JTB| , I7 |, where 
no A dependence is seen until pr exceeds 2 GeV/c. 

On the basis of these arguments we adopt the following position. Firstly, the magnitude 
of k t is small, since heavy quark suffers less deflection; we take it to be of order /i, which 
according to /i 2 = 4iia s T 2 gives /i ~ 0.4 GeV for ot s = 0.3 and T = 200 MeV. Secondly, 
k t is expected to be independent of L, although a gentle increase with L cannot be ruled 
out on firm theoretical ground. Since the recombination of c with a q to form a D after 
passing through the QM does not increase p t beyond the k t gained, we arrive at the result 
that \\ + A 2 , ~ A 2 2 . Consequently, we have A ~ 0.5/p, where p is in units of GeV/c. It is 
clear that if we want A to be small, p should not be small. That is why we have set p > 2 
GeV/c. 

2. Hadron Matter. Consider now the scenario where the created cc pair find themselves in 
a densly packed hadronic medium. The recombination of c with q and c with q take place 
rapidly, and it is the D and D that traverse the HM in opposite directions. As before, the 
combined effect on p\ and P2 can be represented by a D of momentum p going through 
a slab of HM of thickness L. Being a low-momentum hadron (2 < p < 5 GeV/c) the D 
interacts strongly with the mainly pionic medium. The number of multiple collisions in L is 
v = n^acnrL, where n T is the pion density. The energy dependence of is not known, but 
its magnitude (in the few mb range) is definitely much greater than partonic cross section. 
Since A 2 + A 2 , is proportional to u, we therefore expect p t to be larger (compared to the QM 
case) and to increase significantly with L. Herein lies the major difference between the two 
medium effects on £ = pt/p. Since A12 is the same for the two media, the difference does not 
depend on it. 

Putting together the above considerations leads to our suggestion for the signature of 
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QM vs. HM. Measure C(p, £). Plot the £ dependence for fixed p and determine the mean 
£. Examine how £ depends on Et and A. For HM, £ should be large and increase with Et 
and A. If, at high enough s, E T and A, quark matter is created, then £ should drop down 
to a low value and become essentially independent of E T and A. This transition from high 
and increasing £ to low and roughly constant £ is the signature of QM formation. In short, 
quarks are smaller than mesons; their difference should be revealed in the measurement of 

I 

It would also be of interest to study the p dependence of £(p) . If a Dn increase with p, 
then £hm would be constant in p, in contradistinction from £qm(p) which decreases with p. 

Another interesting problem arises when the recombination of c (c) with s (s) is consid- 
ered. Since in QM the hadronization occurs outside the quark phase, the width Cd s d s f° r 
D S D S correlation should not differ significantly from £,£>£>■ However, in HM the formation 
of D s and D s occurs inside the medium, and because of the suppressed D s n scattering due 
to the OZI rule at low energy, we expect £d s d s < £dd- The observation of these differences 
would add to our understanding of what occurs in these systems. 

What is described in this paper is for idealized systems. In reality the system may be 
much more complicated and the transition of the behavior of £ may be very gradual. If 
so, it is unlikely that any other signature would be clear-cut, since it is the system itself 
that is not sufficiently distinctive. Our proposal deals with the nature of the matter probed, 
independent of other inessential complications, such as the absolute normalization of the 
charm production rate, the precise value of the initial temperature, or the validity of pQCD. 
To have a proper theoretical treatment of the problem is essential ultimately, but for now 
the need for a distinctive experimental signature seems to be more urgent. Even if the 
suggested signature turns out to be ineffective because of the complexity not considered 
in this initial investigation, charm correlation should nevertheless reveal much information 
about heavy-ion collisions not available so far. 

I am grateful to X.N. Wang for helpful discussions. This work was supported, in part, 
by the U.S. Department of Energy under grant No. DE-FG06-91ER40637. 
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